Refining DORIS atmospheric drag estimation in preparation of ITRF2008

In preparation of ITRF2008, all geodetic technique services (VLBI, SLR, GPS and DORIS) are generating new solutions based on combination of individual analysis centers solutions. These data reprocessing are based on a selection of models, parameterization and estimation strategy unique to each analysis center and to each technique. While a good agreement can be found for models between groups, thanks to the existence of the IERS conventions, a great diversity still exist for parameter estimation, allowing possible future improvements in this direction. The goal of this study is to focus on the atmospheric drag estimation used to generate the new DORIS/IGN ignwd08 time series prepared for ITRF2008. We develop here a method to inter-compare different processing strategies. In a first step, by analyzing single-satellite solutions for a few weeks of data but for a large number of possible analysis strategies, we demonstrate that estimating drag coefficient more frequently (typically every 1–2 h instead of previously every 4–8 h) for the lowest DORIS satellites (SPOTs and Envisat) provides better geodetic results for station coordinates and polar motion. This new processing strategy also solved earlier problem found when processing DORIS data during intense geomagnetic events, such as geomagnetic storms. Differences between drag estimation strategies can mostly be found during these few specific periods of extreme geomagnetic activity (few days per year). In such a case, when drag coefficient is only estimated every 6 h or less often for single-satellite solution, a significant degradation in station coordinate accuracy can be observed (120 mm vs. 20 mm) and significant biases arose in polar motion estimation (5 mas vs. 0.3 mas). In a second step, we reprocessed a full year of DORIS data (2003) in a standard multi-satellite mode. We were able to provide statistics on a more reliable data set and to strengthen these conclusions. Our proposed DORIS analysis is easy to implement in all software packages and is now already used by several analysis centers of the International DORIS Service (IDS) when submitting reprocessed solutions for ITRF2008.     

基于神经网络模型的襟翼主动控制旋翼减振分析

基于径向基函数(Radial basis function,RBF)神经网络构建了一种带后缘襟翼主动控制(Active controlled flap,ACF)的旋翼振动载荷计算模型。采用正交试验方法确立RBF网络训练样本的输入,在CAMRAD II中计算前飞状态下与训练样本对应的旋翼桨毂六力素,并将主通过频率下的分量作为样本输出,对RBF网络进行离线训练。在此基础上采用多周控制器对被控模型进行振动载荷主动控制。随后以2桨叶4m直径ACF旋翼为例,构建了其桨毂减振分析方法,并对桨毂动载荷各分量的减振效果进行了分析。研究表明,采用正交样本训练的RBF网络能够精确映射襟翼偏角与桨毂振动载荷的非线性关系,施加多周控制后,桨毂垂向振动载荷降低接近50%,其他方向的振动载荷也有不同程度的降低。     

Reconfigurable Flight Control Design for Combat Flying Wing with Multiple Control Surfaces

With control using redundant multiple control surface arrangement and large-deflection drag rudders,a combat flying wing has a higher probability for control surface failures.Therefore,its flight control system must be able to reconfigure after such failures.Considering three types of typical control surface failures(lock-in-place(LIP),loss-of-effectiveness(LOE) and float),flight control reconfiguration characteristic and capability of such aircraft types are analyzed.Because of the control surface redundancy,the aircraft using the dynamic inversion flight control law already has a control allocation block.In this paper,its flight control configuration during the above failures is achieved by modifying this block.It is shown that such a reconfigurable flight control design is valid,through numerical simulations of flight attitude control task.Results indicate that,in the circumstances of control surface failures with limited degree and the degradation of the flying quality level,a combat flying wing adopting this flight control reconfiguration approach based on control allocation could guarantee its flight safety and perform some flight combat missions.     

圆球绕流场的尾涡分析和升阻力研究

采用自编的分块算法程序模拟了雷诺数在20～1000之间的圆球绕流场.在Re=25时捕捉到流动分离,与试验结果一致.对Re=250时稳态非对称尾流区内流体的输运情况进行研究,证实圆球绕流场中纵向对称面的存在.对于非稳态流动,主要研究了阻力和横向力的时间变化规律,并分析了频谱特征.计算发现横向合力主要在横截面上两个很窄的区间内偏移.经过几种不同涡定义方法的优劣测试,本文采用Hunt等人介绍的Q-定义方法对非稳态绕流场的尾涡结构进行了直观描述.     

A Structured Mesh Euler and Interactive Boundary Layer Method for Wing/Body Configurations

To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. For the turbulent boundary-layer, an integral method using Green＇s lag equation is coupled with the outer inviscid flow. A blowing velocity approach is used to simulate the displacement effects of the boundary layer. To predict the aerodynamic drag, it is developed a numerical technique called far-field method that is based on the momentum theorem, in which the total drag is divided into three component drags, i.e. viscous, induced and wave-formed. Consequently, it can provide more physical insight into the drag sources than the often-used surface integral technique. The drag decomposition can be achieved with help of the second law of thermodynamics, which implies that entropy increases and total pressure decreases only across shock wave along a streamline of an inviscid non-isentropic flow. This method has been applied to the DLR-F4 wing/body configuration showing results in good agreement with the wind tunnel data.     

原子钟调相白噪声降噪方法

原子钟的五种噪声反映了它在一定时间间隔内的性能。当取样间隔小于1h，原子钟主要表现为调相白噪声。在地方协调时UTC（NTSC）的实时监控中，不仅要求参考原子时（RTA，用于地方协调时控制的参考）具有很高的长期稳定度，而且要求其短期稳定度亦很高。影响RTA短期稳定度的主要是个台原子钟的调相白噪声。首先讨论了确定原子钟短期噪声的方法，进而利用小波分析方法研究了原子钟凋相白噪声分离与降噪的方法。     

J2项和大气阻力摄动作用下卫星编队构型的演化分析

研究卫星编队构型在J2项摄动和大气阻力摄动作用下的演化进程。从相对运动动力学模型出发,分别加入J2项摄动和大气阻力摄动因素,采用数值积分对所获得的运动方程求解,得到参考星和环绕星在惯性坐标系中的绝对运动规律;将两星的绝对运动规律作差后,利用坐标转换矩阵,将惯性坐标系中两星的相对运动规律转换到Hill坐标系中。由于该方法没有经过任何简化,故研究摄动的影响时不存在误差。最后通过仿真分析,给出了J2项摄动和大气阻力摄动作用下卫星编队构型演化的相关结论,并提出进一步研究的方向。     

基于特征空间变换的运载火箭Pogo模型降阶方法

针对运载火箭Pogo状态空间模型存在维数高和奇异性的问题,开展了Pogo模型的降阶方法研究。基于特征空间变换理论,导出了一般形式的Pogo状态空间模型的解耦形式。该形式与推进系统和结构系统的模态频率相对应,从而通过模态截断或保留感兴趣的模态实现有效的模型降阶。同时,还给出了降阶方法的实数运算公式。在两种不同型号推进系统的火箭Pogo问题中进行了仿真校验。结果表明所提出的降阶方法对奇异和非奇异Pogo状态空间模型都能给出正确的降阶模型,显著地提高计算效率,具有很好的通用性,为Pogo时域仿真和主动抑制提供了合适的模型。     

A salient edges detection algorithm of multi-sensor images and its rapid calculation based on PFCM kernel clustering

Multi-sensor image matching based on salient edges has broad prospect in applications, but it is difficult to extract salient edges of real multi-sensor images with noises fast and accurately by using common algorithms. According to the analysis of the features of salient edges, a novel salient edges detection algorithm and its rapid calculation are proposed based on possibility fuzzy C-means （PFCM） kernel clustering using two-dimensional vectors composed of the values of gray and texture. PFCM clustering can overcome the shortcomings that fuzzy C-means （FCM） cluster- ing is sensitive to noises and possibility C-means （PCM） clustering tends to find identical clusters. On this basis, a method is proposed to improve real-time performance by compressing data sets based on the idea of data reduction in the field of mathematical analysis. In addition, the idea that kernel-space is linearly separable is used to enhance robustness further. Experimental results show that this method extracts salient edges for real multi-sensor images with noises more accurately than the algorithm based on force fields and the FCM algorithm; and the proposed method is on average about 56 times faster than the PFCM algorithm in real time and has better robustness.     

Generalized Kutta–Joukowski theorem for multi-vortex and multi-airfoil flow(a lumped vortex model)

For purpose of easy identification of the role of free vortices on the lift and drag and for purpose of fast or engineering evaluation of forces for each individual body,we will extend in this paper the Kutta–Joukowski(KJ) theorem to the case of inviscid flow with multiple free vortices and multiple airfoils. The major simplification used in this paper is that each airfoil is represented by a lumped vortex,which may hold true when the distances between vortices and bodies are large enough. It is found that the Kutta–Joukowski theorem still holds provided that the local freestream velocity and the circulation of the bound vortex are modified by the induced velocity due to the outside vortices and airfoils. We will demonstrate how to use the present result to identify the role of vortices on the forces according to their position,strength and rotation direction. Moreover,we will apply the present results to a two-cylinder example of Crowdy and the Wagner example to demonstrate how to perform fast force approximation for multi-body and multi-vortex problems. The lumped vortex assumption has the advantage of giving such kinds of approximate results which are very easy to use. The lack of accuracy for such a fast evaluation will be compensated by a rigorous extension,with the lumped vortex assumption removed and with vortex production included,in a forthcoming paper.